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The ratio of the yields of antiprotons to protons in pp collisions has been measured by the 
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and rapidity \y\ < 0.5. The ratio is measured to be R\ y \ <a . 5 = 0.957 ± 0.006(stat.) ± 0.014(syst.) 
at 0.9 TeV and R\ y \ <0 . 5 = 0.991 ± 0.005(stat.) ± 0.0U{syst.) at 7 TeV and it is independent of 
both rapidity and transverse momentum. The results are consistent with the conventional model of 
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baryon-number transport and set stringent limits on any additional contributions to baryon-number 
transfer over very large rapidity intervals in pp collisions. 
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In inelastic non-diffractive proton-proton collisions at 
very high energy, the incoming projectile breaks up into 
several hadrons which emerge after the collision in gen- 
eral under small angles along the original beam direc- 
tion. The deceleration of the incoming proton, or more 
precisely of the conserved baryon number associated with 
the beam particles, is often called "baryon-number trans- 
port" and has been debated theoretically for some time 

One mechanism responsible for baryon-number trans- 
port is the break-up of the proton into a diquark-quark 
configuration The diquark hadronizes after the re- 
action with some longitudinal momentum p z into a new 
particle, which carries the baryon number of the incoming 
proton. This baryon-number transport is usually quan- 
tified in terms of the rapidity loss Ay = ybcam — 2/baryon, 
where f/bcam (j/baryon) is the rapidity of the incoming 
beam (outgoing baryonj^. 

However, diquarks in general retain a large fraction of 
the proton momentum and therefore stay close to beam 
rapidity, typically within one or two units. Therefore, 
additional processes have been proposed to transport 
the baryon number over larger distances in rapidity, in 
particular via purely gluonic exchanges, where the pro- 
ton breaks up into three quarks. The baryon number 
resides with a non-perturbative configuration of gluon 
fields, the so-called "baryon string junction" , which con- 
nects the valence quarks [H, 0] . In this picture, baryon- 
number transport is suppressed exponentially with the 
rapidity interval Ay, proportional to exp [(aj — 1) Ay], 
where aj is identified in the Regge model as the inter- 
cept of the trajectory for the corresponding exchange in 
the t-channel. If the string junction intercept is approxi- 
mated with the one of the standard Reggeon (or meson) , 
aj w 0.5, baryon transport will approach zero with in- 
creasing Ay. If the intercept of the pure stringi unction 
is aj ~ 1, as motivated by perturbative QCD [4], it will 
approach a constant and finite value. 

The LHC, being by far the highest energy proton- 
proton collider, opens the possibility to investigate 
baryon transport over very large rapidity intervals by 
measuring the antiproton-to-proton production ratio at 
midrapidity, R = N p /N p , or equivalently, the proton- 
antiproton asymmetry, A = (N p — Np) / (N p + N p ) . Most 
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" Also at Sezione INFN, Bologna, Italy 

1 The rapidity y is defined as y = 0.5 In [(E + p z ) / (E — p z )]\ ra- 
pidity y = corresponds to longitudinal momentum p z = of 
the baryon in the center-of-mass system and Ay = In (y/s/m p ). 
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of the (anti)protons at midrapidity are created in baryon- 
antibaryon pair production, implying equal yields. Any 
excess of protons over antiprotons is therefore associ- 
ated with the baryon-number transfer from the incoming 
beam. Note that such a study has not been carried out 
in high-energy proton-antiproton colliders (SppS, Teva- 
tron) because of the symmetry of the initial system at 
midrapidity. Model predictions for the ratio R at LHC 
energies range from unity, i.e., no baryon-number trans- 
fer to midrapidity, down to about 0.9 in models where 
the string junction transfer is not suppressed with the 
rapidity interval (aj ~ 1). 

In this letter, we describe the measurement of the 
p/p ratio at midrapidity in non-diffractive pp collisions 
at center-of-mass energies yfs = 0.9 TeV and 7 TeV 
(Ay rj 6.9-8.9), with the ALICE experiment at the LHC. 

ALICE, which is the dedicated heavy-ion detector at 
the LHC, consists of 18 detector sub-systems [8J, l9j. The 
central tracking systems used in the present analysis are 
located inside a solenoidal magnet (B = 0.5 T); they 
are optimized to provide good momentum resolution and 
particle identification (PID) over a broad momentum 
range, up to the highest multiplicities expected for heavy 
ion collisions at the LHC. All detector systems were com- 
missioned and aligned during several months of cosmic- 
ray data-taking in 2008 and 2009 [13, El- 
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Collisions occur inside a beryllium vacuum pipe (3 cm iK 
in radius and 800 fim thick) at the center of the ALICE 
detector. The tracking system in the ALICE central bar- 
rel has full azimuth coverage within the pseudo-rapidity 
window \rj\ < 0.9. The following detector sub-systems 
were used in this analysis: the Inner Tracking System 
(ITS) [ll|, the Time Projection Chamber (TPC) [jj] and'"' 
the VZERO detector H. 

M 109 

The ITS consists of six cylindrical layers of silicon de- 
tectors with radii of 3.9/7.6 cm (Silicon Pixel Detectors-^ . 
SPD), 15.0/23.9 cm (Silicon Drift Detectors-SDD) and m 
38/43 cm (Silicon Strip Detectors-SSD). They provide 
full azimuth coverage for tracks matching the acceptance 
of the TPC (M < 0.9). ^ 

The TPC is the main tracking detector of the central U6 
barrel. The detector is cylindrical in shape with an active 117 
volume of inner radius 85 cm, outer radius of 250 cm and 118 
an overall length along the beam direction of 500 cm. 

Finally, the VZERO detector consists of two arrays of 120 
32 scintillators each, which are placed around the beam m 
pipe on either side of the interaction region) at z = 3.3 m 122 
and z — —0.9 m, covering the pseudorapidity ranges 123 
2.8 < r\ < 5.1 and -3.7 < rj < —1.7, respectively [13]. A 124 
detailed description of the ALICE detectors, its compo- 125 
nents, and their performance can be found in [8]. U6 

Data from 2.8 (y/s = 0.9 TeV) and 4.2 (Vs = 7 TeV)i2 7 
million pp collisions, recorded during the first LHC runsus 
(December 2009, March-April 2010) were used for thisi2<> 
analysis. The events were recorded with both field po-no 
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FIG. 1. (Color online) The measured ionization per unit 
length as a function of particle momentum (both charges) 
in the TPC gas. The curves correspond to expected energy 
loss [3| for different particle types. The inset shows the mea- 
sured ionization for tracks with 0.99 < p < 1.01 GeV/c. The 
lines are Gaussian fits to the data. 



larities for each energy. The trigger required a hit in 
one of the VZERO counters or in the SPD detector, i.e., 
at least one charged particle anywhere in the 8 units of 
pseudorapidity covered by these trigger detectors [l3j |. 
In addition, the trigger required a coincidence between 
the signals from two beam pick-up counters, one on each 
side of the interaction region, indicating the presence of 
passing bunches. 

Beam-induced background was reduced to a negligible 
level (< 0.01%) with the help of the timing information 
from the VZERO counters [13| and by requiring a re- 
constructed primary vertex (calculated from the SPD) 
within ±1 cm perpendicular to and ±10 cm along the 
beam axis. 

Measurements of momentum and particle identifica- 
tion are performed using information from the TPC de- 
tector, which measures the ionization in the TPC gas 
and the particle trajectory with up to 159 space points. 
In order to ensure a good track quality, a minimum of 
80 clusters was required per track in the TPC and at 
least two hits in the ITS of which at least one is in the 
SPD. In order to reduce the contamination from back- 
ground and secondary tracks (e.g. (anti)protons originat- 
ing from weak hyperon decays or secondary interactions 
in the material), a cut was imposed on the distance 
of closest approach (dca) of the track to the primary 
vertex in the xy (transverse) plane, which varied from 
2.65 to 1.8 mm (2.33 to 1.5 mm for the 7 TeV data) 
for the lowest (0.45 < p t < 0.55 GeV/c) and highest 
(0.95 < pt < 1.05 GeV/c) p t bins, respectively. This 
cut corresponds to 5tr of the measured dca resolution for 
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FIG. 2. The distance of closest approach (dca) distributions 
of p and p for the lowest (left plot) and highest (right plot) 176 
transverse momentum bins. The broad background of protons 177 
at low momentum originates from secondary particles created 178 
in the detector material, whereas the tails for both p and pi79 
at high momentum (and for p at low momentum) arise from lao 
weak hyperon decays. 



each momentum bin. 183 
Particles are identified using their specific ionization 184 
(dE/dx) in the TPC gas [Hj]- Figure □ shows the ioniza- 185 
tion (truncated mean) as a function of par ticle momen- 186 
turn together with the expected curves [14[ for different 187 
particle species. The inset shows the measured dE/dx for 188 
tracks in the momentum range 0.99 < p < 1.01 GeV/c 189 
with clearly separated peaks for (anti)protons and lighter 190 
particles. The dE/dx resolution of the TPC is 5, de- 191 
pending slightly on the number of TPC clusters and the 192 
track inclination angle. For this analysis, (anti)protons 193 
were selected within a band of ±3cr around the expected 194 
value. 

In order to assure uniform geometrical acceptance, 196 
high reconstruction efficiency and unambiguous proton 197 
identification, we restrict the analysis to protons and an- 198 
tiprotons in the rapidity range |y| < 0.5 and the momen- 199 
turn range 0.45 < p < 1.05 GeV/c. The contamination 200 
of the proton sample with electrons or pions and kaons is 201 
negligible (< 0.1%) even at the highest momentum bins, 202 
and in addition essentially charge symmetric. 203 

Most instrumental effects associated with the accep- 204 
tance, reconstruction efficiency, and resolution are iden- 205 
tical for primary protons and anti-protons and therefore 206 
cancel in the ratio. However, because of significant dif- 207 
ferences in the relevant cross sections, anti-protons are 208 
more likely than protons to be absorbed or elastically 
scattereco within the detector, and a non negligible back- 



ground in the proton sample arises from secondary inter- 
actions in the beam pipe and inner layers of the detector. 

In order to correct for the difference between p-A and 
p-A elastic and inelastic reactions in the detector mate- 
rial, detailed Monte Carlo simulations based on GEANT3 
[lj| and FLUKA [l6| were performed. These corrections 
rely in particular on the proper description of the interac- 
tion cross sections used as input by the transport models. 
These values were therefore compared with experimental 
measurements 

[H1G3. While p-A cross sections are sim- 
ilar in both models and in agreement with existing data, 
GEANT3 (as well as the current version of GEANT4) 
significantly overestimates the measured inelastic cross 
sections for antiprotons in the relevant momentum range 
by about a factor of two, whereas FLUKA describes the 
data very well. Concerning elastic scattering, where only 
a limited data set is available for comparison, GEANT3 
cross sections are about 25% above FLUKA, the latter 
being again closer to the measurements. We therefore 
used the FLUKA results to account for the difference of 
p and p cross sections, which amount to a correction of 
the p/p ratio by 8% and 3.5% for absorption and elastic 
scattering, respectively. 

The contamination of the proton sample due to sec- 
ondaries originating from interactions with the detector 
material was directly measured with the data and sub- 
tracted. Most of these background tracks do not point 
back to the interaction vertex and can therefore be ex- 
cluded with a dca cut. Figure [5] shows the dca distri- 
butions of p and p for the lowest (left panel) and the 
highest (right panel) transverse momentum bins. Sec- 
ondary protons are clearly visible in the left plot due to 
their wide dca distribution. At higher momenta the back- 
ground of secondary protons becomes very small. The 
remaining tails visible in the dca distributions are due to 
(anti)protons originating from weak decays. The back- 
ground of secondary protons, which remains after the 
dca cut under the peak of primaries, is subtracted by de- 
termining its shape from Monte Carlo simulations and 
adjusting the amount to the data at large values of the 
dca. This correction is calculated and applied differen- 
tially as a function of y and pt] it varies between 14% for 
the lowest and less than 0.3% for the highest transverse 
momentum bins. 

The contamination coming from feed-down (i.e., 
(anti)protons originating from the weak decay of A and 
A) was subtracted in a similar way by parametrization 
and fitting to the data of the respective simulated dca dis- 
tributions. This correction ranges from 20% to 12% for 
the lowest and highest pt bins, respectively. 



2 Particles undergoing elastic scattering in the inner detectors can 



still be reconstructed in the TPC but the corresponding ITS hits 
will in general not be associated to the track if the scattering 
angle is large. 
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TABLE I. Systematic uncertainties of the p/p ratio. 



Systematic Uncertainty 


Material budget 


0.5% 


Absorption cross section 


0.8% 


Elastic cross section 


0.8% 


Analysis cuts 


0.4% 


Corrections (secondaries/feed-down) 


0.6% 


Total 


1.4% 



209 The main sources of systematic uncertainties are the 

210 detector material budget, the (anti)proton reaction cross 
2n section, the subtraction of secondary protons and the ac- 

212 curacy of the detector response simulations (see Table [I]). 

213 The amount of material in the central part of ALICE 

214 is very low, corresponding to about 10% of a radiation 

215 length on average between the vertex and the active vol- 

216 ume of the TPC. It has been studied with collision data 

217 and adjusted in the simulation based on the analysis of 

218 photon conversions. The current simulation reproduces 

219 the amount and spatial distribution of reconstructed con- 

220 version points in great detail, with a relative accuracy of 

221 a few percent. Based on these studies, we assign a sys- 

222 tematic uncertainty of 7% to the material budget. By 

223 changing the material in the simulation by this amount, 

224 we find a variation of the final ratio R of less than 0.5%. 

225 The experimentally measured p-A reaction cross sec- 

226 tions are determined with a typical accuracy better than 

227 5% [iTJ. We assign a 10% uncertainty to the absorption 
22s correction as calculated with FLUKA, which leads to a 

229 0.8% uncertainty in the ratio R. By comparing GEANT3 

230 with FLUKA and with the experimentally measured elas- 

231 tic cross-sections, the corresponding uncertainty was es- 

232 timated to be 0.8%, which corresponds to the difference 

233 between the correction factors calculated with the two 

234 models. 

235 By changing the event selection, analysis cuts and 

236 track quality requirements within reasonable ranges, we 

237 find a maximum deviation of the results of 0.4%, which 

238 we assign as systematic uncertainty to the accuracy of 

239 the detector simulation and analysis corrections. 

240 The uncertainty resulting from the subtraction of sec- 

241 ondary protons and from the feed-down corrections was 1 

242 estimated to be 0.6% by using different functional forms 1 

243 for the background subtraction and for the contribution 257 

244 of the hyperon decay products. 258 

245 The contribution of diffractive reactions to our finals 

246 event sample was studied with different event generators 2 '" 

247 and was found to be less than 3%, resulting into a negligi- 261 

248 ble contribution (< 0.1%) to the systematic uncertainty. 2 <> 2 

249 Finally, the complete analysis was repeated using only 2 " 

250 TPC information (i.e., without using any of the ITS de-264 

251 tectors). The resulting difference was negligible at botIi265 

252 energies (< 0.1%). 266 

253 Table U summarizes the contribution to the system-267 

254 atic uncertainty from all the different sources. The totahes 
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FIG. 3. (Color online) The p t dependence of the p/p ratio in- 
tegrated over \y\ < 0.5 for pp collisions at y/s = 0.9 TeV (top) 
and y/s = 7 TeV (bottom). Only statistical errors are shown 
for the data; the width of the Monte Carlo bands indicates 
the statistical uncertainty of the simulation results. 



systematic uncertainty is identical for both energies and 
amounts to 1.4%. 

The final, feed-down corrected p/p ratio R inte- 
grated within our rapidity and pt acceptance rises from 
%|<o.5 = 0.957 ± 0.006(stat.) ± O.OU(syst) at V« = 
0.9 TeV to %| <0 .5 = 0.991 ± 0.005(s£at.) ± O.OU(syst) 
at y/s = 7 TeV. The difference in the p/p ratio, 0.034 ± 
0.008(stat.), is significant because the systematic errors 
at both energies are fully correlated. 

Within statistical errors, the measured ratio R shows 
no dependence on transverse momentum (Fig. or ra- 
pidity (data not shown) . The ratio is also independent of 
momentum and rapidity for all generators in our accep- 
tance, with the exception of HIJING/B, which predicts 
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FIG. 4. (Color online) Central rapidity p/p ratio as a function 312 
of the rapidity interval Ay (lower axis) and center-of-mass 313 
energy (upper axis). Error bars correspond to the quadratic 
sum of statistical and systematic uncertainties for the RHIC 
and LHC measurements and to statistical errors otherwise. 



269 a decrease with increasing transverse momentum for thesis 

270 lower energy. 319 

271 The data are compared with various model predic-320 

272 tions for pp collisions [H 0, EH in Table [TT] (integrateds2i 

273 values) and Fig. [31 The analytical QGSM model does322 

274 not predict the pt dependence and is therefore not in-323 

275 eluded in Fig. [3] For both energies, two of the PYTHIA324 

276 tunes [lj| (ATLAS-CSC and Perugia-0) as well as the32 5 

277 version of Quark-Gluon String Model (QGSM) with the326 
27s value of the string junction intercept aj = 0.5 Q de-327 

279 scribe the experimental values well, whereas QGSM with- 32s 

280 out string junctions (e = 0, e is a parameter propor-329 

281 tional to the probability of the string-junction exchange)33o 

282 is slightly above the data. HIJING/B [7], unlike thes3i 

283 above models, includes a particular implementation of332 

284 gluonic string junctions to enhance baryon-number trans-333 

285 fer. This model underestimates the experimental results, 334 

286 in particular at the lower LHC energy. Also, QGSM335 

287 with a value of the junction intercept aj = 0.9 [6j pre-336 

288 diets a smaller ratio, as does the Perugia-SOFT tune of 337 

289 PYTHIA, which also includes enhanced baryon transfeo 

290 Figure [4] shows a compilation of central rapidity mea- j39 

291 surements of the ratio R in pp collisions as a function 

1 1 340 

292 of center-of-mass energy (upper axis) and the rapidity 

293 interval Ay (lower axis). The ALICE measurements cor- 342 

294 respond to Ay = 6.87 and Ay = 8.92 for the two energies, 43 

344 
345 



3 We have checked that baryon transfer is the main reason for the 346 
different p/p ratios predicted by the models; the absolute yield 347 
of (anti)protons in our acceptance, which is dominated by pair348 
production, is reproduced by the models to within ±20%. 



whereas the lower energy data points are taken from [20I — 
|22| . The p/p ratio rises from 0.25 and 0.3 at the SPS and 
the lowest ISR energy, respectively, to a value of about 
0.8 at y/s = 200 GeV, indicating that a substantial frac- 
tion of the baryon number associated with the beam par- 
ticles is transported over rapidity intervals of up to five 
units. 

Although our measured midrapidity ratio R at y/s = 
0.9 TeV is close to unity, there is still a small but sig- 
nificant excess of protons over antiprotons correspond- 
ing to a p-p asymmetry of A = 0.022 ± 0.003(stai.) ± 
0.007(syst.). On the other hand, the ratio at y/s = 7 TeV 
is consistent with unity (A = 0.005 ± 0.003(stai.) ± 

0. 007 (syst.)), which sets a stringent limit on the amount 
of baryon transport over 9 units in rapidity. The exis- 
tence of a large value for the asymmetry even at infinite 
energy, which has been predicted to be A — 0.035 using 
aj = 1 [4], is therefore excluded. 

A rough approximation of the Ay dependence of the 
ratio R can be derived in the Regge model, where 
baryon pair production at very high energy is governed 
by Pomeron exchange and baryon transport by string- 
junction exchange [5j. In this case the p/p ratio takes 
the simple form 1/R = 1 + Cexp[(aj — ap)Ay]. We 
have fitted such a function to the data, using as value 
for the Pomeron intercept ap = 1.2 [23j and aj = 0.5, 
whereas C, which determines the relative contributions of 
the two diagrams, is adjusted to the measurements from 
ISR, RHIC, and LHC. The fit, shown in Fig. gj gives 
a reasonable description of the data with only one free 
parameter (C), except at lower energies, where contribu- 
tions of other diagrams cannot be neglected Q ■ Adding a 
second string junction diagram with a larger intercept |4[ , 

1. e., 1/R= l + Cexp[(aj-a P )Ay]+C'exp[(aj/-ap)Ay] 
with a j/ = 1, does not improve the quality of the fit 
and its contribution is compatible with zero (C ~ 10, 
C' « — 0.1 ± 0.1). In a similar spirit, our data could 
also be used to constrain other Regge-model inspired de- 
scriptions of baryon asymmetry, for example when the 
string-junction exchange is replaced by the "odderon" , 
which is the analogue of the Pomeron with odd C-parity; 
see [6j. 

In summary, we have measured the ratio of antipro- 
ton to proton production in the ALICE experiment at 
the CERN LHC collider at yft = 0.9 and y/s = 7 TeV. 
Within our acceptance region (\y\ < 0.5, 0.45 < pt < 
1.05 GeV/c), the ratio of antiproton-to-proton yields 
rises from %|< .5 = 0.957 ± 0.006(stai.) ± 0.014(sysr;.) 
at 0.9 to a value close to unity -Riyi<o.5 = 0.991 ± 
0.005(stat.) ± 0.0U(syst.) at 7 TeV. The p/p ratio is 
independent of both rapidity and transverse momen- 
tum. These results are consistent with standard models 
of baryon-number transport and set tight limits on any 
additional contributions to baryon-number transfer over 
very large rapidity intervals in pp collisions. 
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TABLE II. The measured central rapidity p/p ratio compared to the predictions of different models (the statistical uncertainty 
in the models is less than 0.005). The quoted errors for the ALICE points are the quadratic sum of statistical and systematic 
uncertainties. 



Energy [TcV] 


0.9 


7 


ALICE 


0.957 ±0.015 


0.991 ±0.015 


ATLAS-CSC Tune (306) 


0.96 


1.0 


PYTHIA Perugia-0 Tune (320) 


0.95 


1.0 


Perugia-SOFT Tune (322) 


0.88 


0.94 


e = 


0.98 


1.0 


QGSM e = 0.076, aj = 0.5 


0.96 


0.99 


e = 0.024, aj = 0.9 


0.89 


0.95 


HIJING/B 


0.83 


0.97 
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